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INTRODUCTION 

V  . 

7  Largely  as  a  consequence  of  advances  in  computer  technology 
the  use  of  time  series  analysis  methodology  to  investigate 
correlations  of  brain  electrical  activity  with  clinical  disorders 
has  become  an  area  of  intense  effort  in  many  EEG  research  lab¬ 
oratories.  The  primary  objective  of  many  of  these  investigations 
is  the  detection  of  EEG  time  series  properties  which  offer  a 
quantitative  basis  for  aiding  diagnosis  of  mental  illness  or  brain 
injury.  This  chapter  describes  the  special  eiectrophysiologieal 
tests  which  we  have  been  developing  at  TRIMS*  for  application  to 
research  projects  concerned  with  the  detection  of  EEG  properties 
which  correlate  with  particular  abnormalities  such  as  uncontrolled 
violent  behavior,  learning  disabilities  and  epilepsy. 

SilAIN  SPIKING  ' 

intermittent  deep  bruin  electrical  spiking  as  well  as  scalp  ; 

-  EEG  spiking  has  been  iwi  Seated  ir;  certain  brain  and  behavioral 

1 1 ) 

disorders  such  as  epilepsy  *  ana  uncontrolled  violent  behavior. 

The  correlation  of  abnormal  deep  brain  electrical  spiking  activity 
with  violent  'behavior  has  been  demonstrated  in  nonhuman  primate 
studies by  employing  invasive  methods  wnich  involve  the  surgical 
implantation  of  electrodes  and  subsequent  analysis  of  the  electri-  j 


cal  activity  recorded  from  the  implanted  deep  brain  structures. 
In  light  of  those  results,  a  significant  advancement  in  the 
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diagnosis  of  abnormal  brain  activity,  especially  that  deep  brain 
activity  associated  with  persons  exhibiting  uncontrolled  violent 
behavior,  would  be  achieved  by  the  development  of  detection 
methods  that  are  noninvasive  and  therefore  applicable  in  ordinary 
clinical  BEG  settings. 

Our  HEG  research  on  noninvasive  detection  was  stimulated 
primarily  by  the  initial  finding  of  complex  patterns  of  consistent 
waveshape  in  scalp  EEG  which  were  time-locked  to  spikes  recorded 
from  electrodes  implanted  in  deep  brain  structures  of  rhesus 
monkeys.  These  studies  have  also  shown  that  such  scalp  corre¬ 
lates  of  deep  spiking  can  be  detected  even  in  severe  EEG  noise 
backgrounds  by  the  application  of  digital  filters  appropriately 
designed  to  minimize  the  effects  of  unwanted  EEG  background 
activity,  or  by  special  application  of  cepstral  methods  in  cases 
where  digital  filters  for  pattern  recognition  ore  not  suitable 
because  the  pattern  to  be  detected  is  not  known  a' priori.  The 
analytical  methods  and  their  potential  applications  are  described 
in  references^  ^  ^ . 

DETECTION  METHODS 

The  digital  filtering  procedures  developed  for  detecting 

,-.4. 

scalp  correlates  of  deep  spiking  wer#  based  on  an  analysis  of 
monkey  and  human  scalp  EEG  data  obtained  from  research  projects 
where  simultaneous  recordings  from  deep  brain  structures  were 
available.  Using  the  deop  spike  as  the  trigger  for  averaging 


scalp  EEG  activity,  it  was  observed  that  transient  slow  wave 

activity  frequently  appeared  in  scalp  activity  at  the  same  time 

that  a  spike  occurred  at  depth.  The  waveshape  of  this  transient 

activity  was  usually  distorted  by  the  presence  of  noise; 

therefore  averaging  procedures  were  used  to  achieve  a  better  | 

estimate  of  the  transient  waveshape.  The  power  spectral  density  | 

of  the  scalp  EEG  background  activity  was  also  estimated  in  order  ’ 

to  appropriately  weight  the  spectral  components  in  the  transient 

waveshape  obtained  by  averaging.  The  digital  filter  derived  from 

the  spectral  estimates  of  both  the  transient  waveshape  and  the 

noise  was  employed  as  a  detector  which  looks  at  scalp  activity 

and  reports  on  the  presence  or  nonpresence  of  transient  patterns 

which  match  the  characteristics  of  the  digital  filter,  it  is 

interesting  to  note  that  this  coincidence  of  deep  spiking  and 

transient  EEG  slowing  implies  that  pathological  sharp  spiking  j 

activity  at  depth  produces  slow  wave  activity  at  the  surface.  | 

| 

This  is  consistent  with  clinical  EEC  criteria  which  consider 
focal  alow  activity  to  be  an  abnormal  indication.  , 

The  procedure  for  evaluating  candidate  digital  filters  was  j 
based  on  the  number  of  spikes  detected  in  normal  sub joe to  as  1 

compared  to  the  numbor  of  spikes  detected  in  the  recording  of  ? 

mentally  ill  subjects,  in  the  initial  evaluation  a  comparison 
was  made  of  the  incidence  of  spikes  in  normals'  and  in  violent 
subjects,  based  on  a  Poisson  model  for  random  spiking.  In  this 
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;  model  the  number  of  spikes  detected  over  a  given  length  of  record¬ 
ing  is  compared  with  the  expected  number  derived  frati  data  on 
normals.  The  normal  control  is  used  to  test  the  hypothesis  that 
a  given  record  represents  the  EEG  of  a  normal  subject  under  the 
assumption  that  spikes  in  normal  subjects  are  uniformly  randomly  j 
distributed.  The  methods  of  analysis  underlying  these  evaluation  * 
procedures  depend  on  the  performance  characteristics  of  the  digital 
filter  as  a  detector,  as  well  as  on  the  statistical  model  for 
evaluating  the  significance  of  the  number  of  spikes  detected. 

These  methods  are  described  in  the  next  section. 

COMPUTER  IMPLEMENTATION 

H  By  virtue  of  our  con$)Uter  configuration,  Fourier  serios 
f methods  (rather  than  matrix  inversion  methods)  wore  used  to  design 
as  well  as  evaluate  the  digital  filter.  The  Fourier  methods 
for  the  design  of  the  optimum  filter  proceed  as  follows t 
I  1.  Obtain  (a)  the  complex  Fourier  series  of  the  candidate 
|  transient  patterns  and,  (b)  the  power  spectral  density  of 
the  signal  plus  noise;  that  is,  the  background  EEG  signal. 

2.  Divide  the  conjugate  complex  Fourier  series  of  the  j 

transient  pattern  by  the  power  spectral  density  of  tho  j 

signal  plus  noisu. 

3.  Take  the  inverse  transform  of  the  Fourier  series  obtained 
in  Step  2  which  gives  u  discretely  sampled  time  function 


that  is  the  desired  template  or  matched  filter. 


In  addition  to  performing  the  above  operations,  the  computer 
in  our  laboratory  is  also  capable  of  performing  running  convolu¬ 
tion.  This  allows  continuous  digital  filtering  of  the  scalp  EEG 
to  rapidly  evaluate  the  performance  characteristics  of  a  candidate 
digital  filter  as  a  detector  of  abnormal  transient  activity. 

The  above  analytical  procedures  refer  to  the  detection  of 
spiJce  induced  events,  but  it  is  necessary  to  assign  sane  signifi¬ 
cance  to  the  number  of  events  detected  in  terms  of  background 
activity  and  artifacts  that  produce  false  spike  indications.  The 
major  difficulty  which  presents  itself  in  physiological  signal 
studies  of  this  type  arises  from  the  fact  that  any  given  signal 
characteristic  such  as  a  spike  can  and  usually  does  appear  due  to 
random  background  effects.  The  problem  thon  becomes  one  of 

determining  whether  the  appearance  of  this  signal  characteristic 

* 

is  due  to  background  effects  or  to  some  inheront  neurophysiolog¬ 
ical  abnormality  in  the  EEG  being  analyzed.  The  analysis  and 
evaluation  rationale  are  straightforward  if  it  is  assumed  that 
the  time  interval  between  spike  indications  is  random  and  uni¬ 
formly  distributed  duo  to  background  activity  in  the  normal  EEG. 
With  the  foregoing  assumptions,  the  analysis  proceeds  as  follows: 
let  p  «  probability  of  oxactly  n  spikes  occurring 

in  time  t  due  to  EEG  background  activity  in 
normal  subject:)  and, 


p  (t+At)  =  probability  of  exactly  n  spikes  occurring 

in  time  t+kt  due  to  EEG  background  activity 
in  normal  subjects  and, 

X  =•  average  rate  at  which  spikes  occur  in  the  EEGs 
of  normal  control  subjects; 

then  p  (t+At)  =*  p(t)  (l-\At)  +  r>  .(tjAAt.  (1) 

Equation  (1}  states  (a)  that  the  probability  of  exactly  n 
spikes  occurring  over  t+At  is  equal  to  the  probability  that  n 
spikes  occur  over  time  t  and  none  in  At,  plus  the  probability  that 
exactly  (n-1)  spikes  occur  over  time  t  and  exactly  one  in  At  and 
that  the  time  increment  At  is  so  small  that  (b)  the  probability  of 
two  or  more  spikes  occurring  during  At  is  zero,  and  (c)  that  the 
probability  of  one  spike  occurring  at  At  is  XAt,  Rearranging  terms 
in  Equation  (1) ,  and  passing  to  the  limit  as  At  •+  0  gives: 
dp  (t) 

-Sr  *  Xpn(t)  “  Xpn-l(t)-  i2) 

The  solution  of  this  difference-differential  equation  is  the 
Poisson  probability  density: 


l’nn> 


at)"  -xt 

n\ 


(3) 


The  average  br  expected  number  of  spike  occurrences,  n,  during 
time  t  is  given  bys 
H  *  Xt. 

The  variance  of  spike  occurrences  is  also  equal  to  \t,  sot 
0*  -  rt. 
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If  the  actual  number  of  spikes  detected  in  an  EEG  record  of 
length  t  is  N,  then  to  test  significance  we  need  to  determine  the 
likelihood  that  N  or  more  spikes  could  occur  in  a  normal  EEG 
record.  This  is  given  by  the  cumulative  distribution: 


P  (t)  =  £  fXtJ 

7  Nf  (\t)n  -Ai 
*  2  -  2  * — — r—  e 


»»0 


nl 


(4) 


If  #  is  large  compared  to  n(s\t)3  then  the  likelihood  computed 
from  Equation  (4)  is  small.  If  this  likelihood  is  sufficiently 
small,  then  we  reject  the  hypothesis  that  an  EEG  record  containing 
N  spikes  in  a  time  t  is  a  normal  record. 

For  several  values  of  spike-count  expectation,  tho  level  of 
significance  associated  with  this  hypothesis  for  a  record  with  »V 
spike  indications  can  be  obtained  from  plots  of  the  cumulative 
distribution  shown  in  Figure  1.  For  example,  if  tho  expected 
nurabor  of  spikes  over  a  given  length  of  a  normal  record  is  6,  then  | 
the  plot  shows  that  the  probability  of  18  spikes  occurring  in  a 
normal  record  is  0.0001.  Therefore,  the  hypothesis  that  an  EEG  in 
which  18  spikes  are  detected  represents  a  normal  subject  is  rejected 

l 

at  the  0.0001  level.  Tho  detection  of  12  spikes  would  allow  re-  j 
jection  of  the  hypothesis  at  the  0.02  level. 

It  should  be  pointed  out  that  such  statistical  modeling  of 
multiple  detections  over  long  EEC  records  Is  essential  because  »*f 
the  false  signals,  artifacts,  and  the  many  uncontrollable  sources 
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of  physiological  interference  that  plague  the  evaluation  of  EEGs. 
These  statistical  procedures  follow  and  supplement  the  digital 
filtering  procedures  used  in  designing  the  detector,  as  outlined  in 
the  previous  sections  and  described  in  our  publications. 
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DETECTION  OF  SPIKE  CORRELATED  SCALP  TRANSIENTS  OF  UNKNOWN  WAVESHAPE 
An  alternate  approach  to  the  deep  spike  noninvasive  detection 
problem  is  required  when  the  recurring  scalp  EEG 'transient  pattern 
is  of  unknown  waveshape.  Since  averaging  methods  for  visualizing  i 
a  recurrent  waveshape  in  noisy  EEG  require  that  the  averaging  pro-  J 

t 

cess  be  synchronized  by  the  deep  spike  event  which  can  be  detected 

< 

only  by  invasive  methods,  the  waveshape  of  the  recurrent  transient 
is  frequently  unknown.  Under  these  conditions,  the  methodology^  ^ 
for  detecting  the  presence  of  a  recurrent  complex  transient  in 
scalp  recorded  brain  electrical  activity  is  based  on  the  application 
of  deconvolution  procedures  aa  described  below, 

A  recurrent  transient  waveform  in  the  KEG  can  be  represented 
as  follows* 

»  , 

AY t)  -  >;  AltXU- xy)  +  Nit)  1  (A) 

k*0  K  * 

Alternatively,  (S)  may  be  written  as  a  convolution  product: 

M 

Sit)  •  A  it)  *  t  A.  6(t~tJ  f  Nit)  (6) 

k*d 

where  *  designates  convolution,  and 

X(t)  “  intermittent  pattern;  i.o.,  waveshape  of  recurrent 
transient 

Nit)  **  EEG  background  activity 
-  6(t^X^)  *  Dirac  delta  function  at  x ^ 

The  deconvolution  of  the  convolution  factors  in  equation  (6)  is 
accomplished  in  several  steps.  First,  the  Fourier  tr.uiM •*! 


» 


(6)  gives  the  algebraic  product  of  the  individual  Fourier  trans-  i 
forms  of  the  intermittent  pattern  and  the  set  of  delta  functions.  | 
This  suggests  the  use  of  cepstral  analysis  which  involves  com¬ 
putation  of  the  logarithm  of  the  Fourier  transform  as  a  second 
step  and,  as  a  third  step,  computation  of  the  inverse  Fourier  | 

transform  of  this  result  to  produce  a  function  called  the  cep-  1 

strum.  The  properties  of  the  cepstrum  will  reveal  the  presence 
of  a  recurrent  pattern  in  the  EEG  by  virtue  of  spikes  which  will 
appear  in  the  cepstrum  when  two  or  more  recurrences  of  the  pattern 
are  embedded  in  the  EEG  epoch  analyzed. 

If  the  wavefora  characteristics  are  of  interest,  then  this 
methodology  can  also  be  used  to  determine  the  shape  of  the  tran¬ 
sient  pattern.  This  is  accomplished  by  smoothing  the  cepstrum  to 
eliminate  the  spikes,  and  then  reversing  ell  the  transformations 
used  to  produce  the  cepstrum.  However,  this  is  a  difficult  eom- 

i 

putatioual  problem  and  it  is  possible  to  circumvent  these  procedures 
if  one  is  not  interested  in  ascertaining  the  shape  of  the  pattern, 
knit  simply  in  detecting  whether  a  recurrent  pattern  is  contained  , 
within  the  time  epoch  analyzed.  If  at  least  two  patterns  are 
captured  in  the  data  epoch,  then  analysis  shows  that  the  (tower  ji 
spectral  density  (PSD)  will  contain  ripples  which  are  attributable  * 
to  the  presence  of  the  recurring  pattern.  The  assumption  under¬ 
lying  the  utility  of  this  approach  is  that  the  background  REG,  in 
the  absence  of  a  recurrent  transient  pattern,  will  possess  a 
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smooth  or  unrippled  PSD.  Figures  2  and  3  demonstrate  that  this 
assumption  holds  for  the  EEG  data  recorded  under  "eyes  open" 
conditions  from  occipital  leads  during  an  experiment  in  which 
transients  were  introduced  into  the  background  EEG  by  intermittent 
visual  stimulation.  The  figures  show  that  the  PSD  for  the  no 
stimulus  condition  is  smooth  (Figure  2)  while  the  PSD  for  the 
stimulus  condition  exhibits  ripples  (Figure  3)  whose  peaks  arc 
separated  by  the  reciprocal  of  the  stimulus  intorval. 

In  summary,  the  above  results  demonstrate)  that  PSD  analysis 
of  sufficient  frequency  resolution  to  resolve  ripples  may  provide 
a  tool  for  noninvasively  diagnosing  illnesses  in  which  deep 
brain  electrical  spiking  may  be  a  factor.  More  generally,  the 
analytical  methods  described  in  this  section  provide  a  basis  for 
investigating  the  clinical  implications  of  weak  recurrent  transient! 
which  are  embedded  in  EEC  background  and  therefore  usually  not 
discernible  by  visual  inspection  of  the  EEC*  time  series. 


EEG  SPECTRAL  COHERENCE 


Spectral  coherence  analysis  of  the  EEG  provides  a  frequency  1 
dependent  measure  of  shared*  electrophysiological  activity.  Thus,  •, 
if  linearly  related  (i.e.,  coherent)  electrophysiological  activity 
between  two  EEG  channels  is  present  in  a  restricted  portion  of  the 

! 

frequency  spectrum  while  the  remainder  of  the  spectrum  contains  | 
activity  which  is  linearly  independent  (i.e.,  incoherent)  then  the 
spectral  coherence  function  by  virtue  of  its  frequency  dependence, 
can  detect  coherent  activity  even  in  situations  where  intense 
levels  of  incoherent  activity  dominate  the  energy  spectrum.  Since 
differences  in  shared  EEG  activity  may  reflect  differences  in 
neural  connectivity  (i.e.,  communication  between  brain  regions) 

(7) 

this  measure  (coherence)  has  been  adopted  by  several  investigators 
as  a  logical  approach  to  the  study  of  brain  function  in  projects 
dealing  with  EEG  correlates  of  cognition  and  learning  disability. 

j 

PROBABILITY  DISTRIBUTION  OF  COHERENCE  ESTIMATES  i 

l 

i 

The  coherence  function  (i.e..  Spectral  Coherence)  is  defined  in 


terms  of  the  normalized  cross-spectrum  of  two  time  series.  The 

. 

cross-spectrum  is  defined  as  the  Fourier  Transform  of  the  cross- 
correlation  function,  viz:  , 


*In  this  context  shared  electrophysiological  activity  among  brain 


regions  is  defined  as  that  activity  *at  a  recording  site  that  is 
related  to  the  activity  at  another  recording. site  through  a  linear 
transformation . 


Denote,  S\(t),  Sift)  =  different  time  series 

<t>nfx)  -  cross-correlation  function  of 

S\(t),  Sift),  where  T  is  the  variable 
time  shift  between  S\  and  S 2,  and  fynfx) 
is  defined  by  the  integral  equation  (7) 


) 


«  —  /J  Sift)  Sift  +  x)dt  (7)  i 

Then  denote,  Fuff)  ~  cross-spectrum  of  S\ft),  Sift)  where  Pn(J') 


is  defined  by  (8) ,  using  the  exponential  form  of  the  Fourier  trans¬ 
formation  . 


Pi  2  (TV  =  (0) 

(where  f  is  frequency  in  hertz  and  T  is  time  shift  in  seconds) 

By  substituting  equation  (7)  into  equation  (8)  and  appropriately 
factoring  the  resulting  double  integral  it  can  be  shown  that 

Pliff)  =  Siff)  Sfz(f)  (c  denotes  complex  conjugate)  (9) 

where:  ~Si(f)  is  the  Fourier  transform  of  S\(t)  and  IPzff)  is  the 
conjugate  Fourier  transform  of  Sift).  Note  that  Siff)  and 

diff)  are  complex  numbers  in  general  and  therefore  so  is  Puff). 
We  may  represent  them  in  polar  form,  as  follows: 

Si(f)  »  riel9‘  l10> 

'H(f)  =  r2e'{0!  (11) 

Substituting  (10)  and  (11)  into  (9)  qives  the  cross-spectral 
density  in  polar  form: 


Pi  2  (f)  =  v\v<ie 


-  vwze 


{where  8  =  0i  -  Q-z) 

While  not  explicitly  shown,  note  that  the  r's  and  8’s  are 
functions  of  j.  Now  to  arrive  at  the  spectral  coherence  function 
we  normalize  equation  (12)  by  dividing  by  rii>2  and  then  averaging 
over  a  number  of  frequencies  and/or  averaging  over  a  number  of  time 
epochs  of  the  two  time  series  being  analyzed. 

Thus  if  the  average  is  taken  over  2N  +  1  discrete  frequencies 
we  may  write  the  spectral  coherence  as 


Q(  f  )  .  - _ i_  y 

CfV  2N+1  e  K 


where  fo  is  the  center  of  the  spectral  window  and  uniform  weighting 
is  used  in  the  window.  If  the  average  is  taken  over  M  time  epochs 
(ensemble  averaging)  we  may  write  the  spectral  coherence  as 

/  :hWk)  -  '4  Jj  (14)  \ 

i 

If  a  combination  of  frequency  and  ensemble  averaging  is  used  then  ] 
spectral  coherence  may  be  written 

1151  , 

Age  these  expressions  assume  the  use  of  uniform  weighting  I 
in  the  spectral  window. 

’  It  should  be  noted  that  for  linearly  independent  signals  the 

■i  ..  ’ 

phase  difference  0C-G4-02i  is  a  random  function  over  both 
JregtfOncy  .and  oasomblo  and  therefor  the  oxpoctod  value  of  coherence  ■ 
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averaged  over  frequency  and/or  time  (ensemble  averaging)  is  zero. 

Thus,  when  coherence  deviates  significantly  from  zero,  one 

de  that  the  two  signal  processes  in  question  are  related 
linear  transformation  over  the  spectral  region  where 
such  significant  deviation  occurs. 

The  Statistical  measure  of  significance  of  a  particular 

§ 

coherence  estimate  depends  on  the  number  of  independent  samples  of 
coherence  which  are  used  in  arriving  at  the  coherence  ostim-'tv:,  and 
upon  the  probability  distribution  function  of  these  coherence 
estimates  as  described  below. 

A  coherence  estimate  is  the  average  computed  from  independent 
samples  rf.  the  normalized  cross-spectrum.  These  samples  can  be 
represented  appoints  on  the  unit  circle  in  the  complex  plane, 
as  illustrated  in  Figure  4.  If  those  points  are  uniformly 

distributed  over  the  unit  circle,  it  is  clear  that  the  expected 
value  of  coherence  lies  at  the  origin  (x  a  y  a  Q),  Thu  normalized 
sampled  or  cross-spoctxcra  may  be  written! 

-i'Q 

xk+^lJk  ~  e  k  *  ooefi^  +  i  oin  0^ 

The  distribution  shown  in  Figure  4  is  equivalent  to  statlny 
that  tho 'relative  phase  values  botwaon  Si  and  S2  at  frequency  / 
are  unifor$lyvr<«idonviy  distributed  over  an  ensemble  of  time  epochs 
and/or  over  a  number  of  frequencies  in  a  spectra)  window  centered 
at  frequency  /  .  The  complex  value  of  coherence  may  bo  written: 


may  concli 


through  a|f 


v 

i 


(16) 


2  *  ,  N 

Coh(f0}  “  n  kl2xk+ 11  n 


Since  x ^  and  y ^  lie  on  the  unit  circle 


**  *  oo8$k 


Vk  = 


the  PDF  (Probability  Distribution  Function)  of  both  x  and  y 
(given  that  0  is  a  uniformly  distributed  random  variable)  have  the 
same  form  and  are  given  by: 

DM  *  ,  D(y)  *  -■■■- — "  (H) 

-a?2  'ir/j  -  y* 

The  corresponding  means,  and  variances  (o\  of)  are  liven  by 

x  =  y  *»  0  (18) 

o*  «  a2,  o  (19) 

H  cVr 

whore  N  *  number  of  samples. 

From  (18)  and  (19)  it  follows  that  the  mean  and  variance  of  the 
coherence  magnitude  C\ooh(f  ) |  S  v)  is  given  by: 

n  =>  0  rjtw 


j  «  —a 

*»  n 


Thus,  to  test  the  hypothesis  that  |t  (t)  and  are  linearlyl 

independent  time  series,  one  examines  the  probability  that  the 
empirically-obtained  moan,  computed  from  N  samples,  deviates  from 
tho  expected  value  (aero  in  this  case) The  number  of  standard 


4m 

< 

deviations  by  which  the  empirical  value  exceeds  the  expected 
value  determines  the  level  of  significance. 

For  example  when  the  coherence  magnitude,  v,  is  obtained  by 
a  combination  of  ensemble  averaging  over  20  epochs  and  frequency 
averaging  over  5  spectral 'components  then 

n  «  so  x  a  *  ioo .  ' 

and  the  resulting  standard  deviation  of  t 

0  »  —  »  — - —  -0.1 

r  sr  rm 

Thus  an  empirically  obtained  value  of  coherence  magnitude 
which  exceeds  0.2  would  be  more  than  two  standard  deviations  from 
the  expected  value  for  independent  signals.  Therefor  it  would  be 
statistically  reasonable  to  conclude  that  the  signals  in  question 
are  not  independent,  but  ruthor  that  they  are  linearly  dependent 
or  coherunt  over  tho  spectral  region  whore  the  coherence  magnitude 

t 

| 

exceeds  0,2.  At  TRIMS  this  method  of  ootiorettce  analysis  has  been  ( 

i 

applied  to  a  pilot  study  of  reading  disabled  children  and  normal 

*  i 

controls .  Our  initial  findings  suggest  that  bilateral  EEG  coherence 
at  frequencies  above  20hs,l*  significantly  lower  in  reading  disabled 
children  than  in  normal  readers,  which  may  be  attributable  to  | 
reduced  sharing  or  communication  between  hemispheres  at  those  ,  (I 
frequencies.  In  any  event  this  EGG  measure  provides  a  basis  for 
comparing  activity  which  visual  examination  of  the  EGG  is  incapablo 
of  discerning  because  tho  lower  frequency  energy  which  tluminufen 


ho  estimate  is 


both  channels  conceals  the  coherent  relationship  at  high  frequen¬ 
cies.  The  application  of  coherence  analysis  to  problems  in 
electrophysiology  is  likely  to  grow  in  importance  because  of  the 
need  to  examine  the  relationships  among  multiple  channel  activity 
which  may  reveal  abnormalities  that  cannot  be  found  in  analysis 
of  single  channel  properties  of  the  EEG. 
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Resting-*  Eyes  Open,  No  Stimulation 


Figure  2.  —  Resting  —  eyes  open  EBG  without  stimulation,  and  corresponding 
Power  Spectral  Density. 


Figure  4.  A  coherence  sample  as  a  point  on  the  unit  circle 
in  the  complex  pUno.  A  coherence  estimate  is  the  average  of  such 
samples  over  a  specified  number  of  frequencies  and/or  a  specified 
number  of  t  ime  ettochn. 


